Artificial electric field in Fermi Liquids 
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Based on the Keldysh formalism, we derive an effective Boltzmann equation for a quasi-particle as- 
sociated with a particular Fermi surface in an interacting Fermi liquid. This provides a many-body 
derivation of Berry curvatures in electron dynamics with spin-orbit coupling, which has received 
much attention in recent years in non-interacting models. As is well-known, the Berry curvature 
in momentum space modifies naive band dynamics via an "artificial magnetic field" in momentum 
space. Our Fermi liquid formulation completes the reinvention of modified band dynamics by intro- 
ducing in addition an artificial electric field, related to Berry curvature in frequency and momentum 
space. We show explicitly how the artificial electric field affects the renormalization factor and trans- 
verse conductivity of interacting (7(1) Fermi liquids with non-degenerate bands. Accordingly, we 
also propose a method of momentum resolved Berry's curvature detection in terms of angle resolved 
photoemission spectroscopy (ARPES). 

PACS numbers: 72.10.Bg, 72.15.Gd, 79.60.-i, 71.10. Ay, 71.10.-w 



Recent experimental developments in spintronics, no- 
tably several measurements of evidently intrinsic anoma- 
lous Hall effect in ferromagnetic metals and observa- 
tions of the spin Hall effect in semiconductors [2, 3], have 
highlighted the importance of intrinsic spin orbit effects 
in conductors. Theoretically, there has been a renewal of 
interest in the fundamental dynamics of Bloch electrons 
in periodic solids taking into account spin orbit inter- 
actions, following ideas pioneered by Karplus and Lut- 
tinger [J], and elegantly formulated by Sundaram and 
Niu [Bj]. In particular, in a non- interacting band pic- 
ture, the low energy Hilbert space is restricted to states 
in the conduction band containing the Fermi surface 
at each quasimomentum k along this surface. This lo- 
cal constraint, like other local constraints in physics Q, 
generates a gauge symmetry in momentum space, cor- 
responding to the freedom to redefine Bloch functions 
by a momentum-dependent phase factor (or SU(2) ma- 
trix in the case of two- fold degenerate bands). Conse- 
quently, a wavepacket analysis of non-interacting Bloch 
states shows that this structure induces artificial mag- 
netic field (a.m.f.), which contributes a Lor entz force in 
k-space, mathematically analogous to the action of the 
real magnetic field b in real i?-space [5j. The a.m.f. 
B a and its associated magnetic gauge field A a are de- 
fined from the periodic part of the Bloch wavefunction 
of conduction bands \u a ): B a j = ieijmdkjA a ,m and 
A a .i = (u^dkiUa)- The fundamental framework of 
Fermi liquid theory suggests that these conclusions per- 
sist in realistic metals, in which the Coulomb interac- 
tions between electrons are not weak. Haldane has ar- 
gued that this is indeed the case for quasiparticles at 
the Fermi surface, with \u a ) replaced by the eigenvector 
of the inverse of the Hermitian part of the retarded ( ad- 
vanced) Green function [8j|. In this letter, we identify an 
additional manifestation of this gauge structure unique 



to interacting Fermi liquids. In particular, because the 
renormalized eigenstates thus introduced are in general 
dependent on frequency u>, it is natural to introduce an 
artiheial electric field (a.e.f.) and electrostatic field as 
£<x,i = i{dujA ati - d ki Aa :0 ) and A afi = (ua\da,u a ) re- 
spectively. These quantities vanish in the absence of in- 
teractions, since the w-dependence of the quasiparticle 
states arises from that of the (collisional) self-energy. 

In the following, we present explicit calculations that 
demonstrate this a.e.f. indeed has physical effects. In 
addition, our Keldysh formulation provides a many-body 
derivation of an effective kinetic equation embodying the 
entire emergent gauge structure, and a systematic pro- 
cedure for actually carrying out the projection into the 
low energy band. We concentrate on metallic materi- 
als lacking either inversion or time-reversal symmetry, in 
which case the bands are non-degenerate and the gauge 
structure is U(l). For such a U(l) Fermi liquid (FL), 
we demonstrate that both a.m.f. and a.e.f., which are 
estimated "on-sheli" (see Eq. (fl2l in combination with 
Eqs. (|9I7[) for their more precise definitions), in fact en- 
ter into the effective equation of motion (EOM) for a 
quasi-particle; 
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£ a thus provides another source of the Lorentz force in 
k-space, in addition to the contribution of B a field identi- 
fied in Ref. [f| It therefore modifies the transverse current 
carried by a quasiparticle, and hence the intrinsic anoma- 
lous Hall conductivity of a conducting ferromagnet to 
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where the sum a is over bands, and n a is the distribution 
function. In particular, Eq. Q can be also proven to be 
consistent with the many-body formulation of the trans- 
verse conductivity more directly by using Ward 
identites which relate the current vertex correction to sin- 
gle quasi-particle Green functions [l2j • As another phys- 
ical consequence of these artificial electromagnetic fields 
(a.e/m.f), we will also show that the linear response of 
the momentum-resolved single-particle density of states 
to applied electromagnetic fields is characterized by these 
artificial fields. Specifically, in the presence of e and b, 
the quasiparticle renormalization factor of a U(l) FL ac- 
quires the following topological corrections: 

Z' = 1 - \ {B a ■ b + E a ■ e) . (3) 

This in principle makes it possible to detect the distribu- 
tion of these "artificial" fields in the momentum space, 
using ARPES experiments. 

We now turn to the derivation of the above results and 
the gauge-invariant effective kinetic equation. We begin 
with a quite general semi-microscopic model in which the 
electronic spectrum is described by a k -p type expansion 
about some (arbitrary) point in the Brillouin zone. One 
may keep as many bands as are deemed close enough 
in energy to be relevant to the physics, and our argu- 
ments do not depend upon the order in the expansion in 
k. Familiar examples would be the multi-band Luttinger 
models for commonly studied semiconductors, in which 
the natural expansion is about the T point. In principle 
one should choose the expansion point so that there are 
no low-energy band crossings intervening between it and 
the Fermi surface, and there could be difficulties in treat- 
ing open Fermi surfaces with the periodic structure of 
the quasimomentum. However, while a fully microscopic 
treatment would be desirable, our approach is neverthe- 
less quite general and moreover we speculate that our fi- 
nal results naturally extend and apply to arbitrary U(l) 
FLs. 

The advantage of this formulation is that we can 
Fourier transform in the usual way to a continu- 
ous real space coordinate r, i.e. k — > — iV r . The 
non-interacting Hamiltonian is thereby expressed 
in terms of the fermionic "envelope fields" ip a {r)'- 
Wo = J2 a , a > I dripl(r)[H (-iV r ,r)] aa ,ip a ,(r), where 
the band index runs from 1 to iVj,. A specific in- 
formation about lattice, orbital and spin-orbit couplings 



are encoded into the matrix structure of [Ho]- We 
will employ a natural form of the interaction, Hi — 

though the form of our results does not depend in detail 
upon this. 

Bearing in mind this Hamiltonian, H = Hq + Hi, we 
can proceed with the Keldysh formalism, beginning with 
an analysis of the spectral function [A(r, t\r' , t')]( n \ n ,\ = 
({Mr,t),^l(r',t')} + ). It obeys 0,0 

[Gq 1 - E HF - a, A] 01 _ ee [L, A]®,_ = [f , ReG] 8 ,_ - 6.(4) 

Here the product (g> denotes the convolution w.r.t. time 
t, space r and band index a, (B®C)(l\l') = J 
(7(1 1 1'). The commutator here is defined by using this 
convolution: [B,6]®- = B ® C - C ® B. The bare 
Green function is as usual defined by Gq 1 (1 1 1') = [id H l- 
^oHV ri) n)] QiQ ^(ii - t[)S(n - r[). E HF (1|T) and 
<t(1|1') are the Hartree-Fock part (temporally instan- 
taneous) and the hermitian part of the collisional self- 
energy (non-instantaneous), while r(l|l') is the anti- 
hermitian part of the collisional self-energy. In a FL, 
the decay rate represented by the [T, ReG]®,- in Eq. (j4|) 
is expected to be 0((uj — fi) 2 ,T 2 ), due to the phase-space 
constraints inherent to FL theory. Because we focus on 
low temperature physics (and indeed require the self- 
energy only to 0(cj — /z), the order that determines the 
non- vanishing electric field at the Fermi surface [HI]), we 
henceforth ignore this lifetime effect, replacing the right 
hand side in Eq. (|U) as indicated. 

Even this dissipationless Keldysh equation is still non- 
trivial and involved, due to the presence of the band in- 
dex, whose effect is the central issue of this work. In the 
following, we will argue a general method of projecting 
out the irrelevant band indices and how to obtain the 
reduced Keldysh (kinetic) equation only for the relevant 
bands. 

The first step is to apply the Wigner transforma- 
tion, so that frequency and momentum (to,q) are intro- 
duced by the Fourier transformation for the relative co- 
ordinate in space and time (t,r) — (t\ — ty,r\ — ry). 
Namely, B(uj,q;T,R) = J drdtB(r u ti\n>, t v ) e -^+^\ 
where the center of mass coordinate in space and time 
X = (T,R) EE { tl+ 2 1 ' » ri+ 2 1 ' ) parameterize phase space 
in combination with Q = (cu,q). Then, following Ref. @, 
the Keldysh equation in terms of Q,X becomes 



- [ L ' k ]-= l 2 [dx ] L,d Q] A] + - i ([d Xj B Xh L,d Qj d Qh A]_- {X k <- Q k }) - {X^Q,}. (5) 

I 



The j and fc-summation over to d are implicit from now on (d is the spatial dimension of our system). The 
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(anti)commutator here is defined in the sense of the prod- 
uct only with respect to band indices. The non-local 
correlation effect which was previously encoded in the 
space-time convolution is now partially taken into ac- 
count via the so-called gradient expansion in 8q8x = 
~d u &T + d Qj . The expansion is justified close to equi- 
librium, since all quantities are independent of X in equi- 
librium. We carry out the gradient expansion to 2nd or- 
der, at which E , B appear in the kinetic equation (see, 
for example, Eq. (fH)) '). 

The Wigner-transformed Lagrangian L, when diago- 
nalized, L d = ■ L ■ U , defines a renormalized energy 
dispersion by its eigenvalue and renormalized Bloch func- 
tions by its eigenbasis. Namely, in this new basis, Eq. ([5]) 
at equilibrium can be satified by those spectral functions 
having no off-diagonal components. The diaognal ele- 
ments of A = U^-A-U are verified a posteriori to become 
a delta-function whose argument is the eigenvalue of L 

A a =A aa =S(L dta ) = S((j-E a (uj ) q)), (6) 

where E a (uj, q) is the a-th eigenenergy of Ho + E HF + a. 

Out of equilibrium, the off-diagonal element also ac- 
quires small but finite weight, which we will take into 
account in a systematic way below. In this new basis, 
the simple derivatives in Eq. ((5]) is replaced by covariant 
ones, 

C/t • (d Qj L) ■ U = d Q l d + A Qj -L d -L d - A Qj , 
C/t • (d Xj A) -U = d Xj A + A Xj ■ A - A ■ A X] , 
Tfl-[L,k]_-U = [L d ,A]_, A Qj = frd Qj U, (7) 

Nb times Nb matrices Aq 6 and A Xj are the renormal- 
ized gauge field, through which we will define a.e/m. f. 
later (see Eqs. (| 12I9() . As the off-diagonal components 
of the spectral function are expected to be small in this 
new basis, they may be eliminated in favor of the diag- 
onal components. Specifically, this occurs because the 
Keldysh equation for the A a p (a ^ (3) contains the term 
^■apAap, with the direct band gap A a p — E a — Ep as 
the sole non-vanishing term at th order in the gradient 
expansion - forcing A a p to vanish in equilibrium. Ac- 
cordingly, A a p is of at most 1st order in the gradient ex- 
pansion, and it can be solved for in terms of the diagonal 
elements A a by iteration. By contrast, the diagonal com- 
ponents of the Keldysh equation appear all at first and 
higher order in the gradient expansion. Substituting this 
iterative solution for A a p into the diagonal equations, we 
can obtain a set of Nb equations for the diagonal parts 
A a = A aa alone, to the desired (2 nd ) order in gradients. 

Finally, we focus on a specific band "a" which con- 
tains a Fermi surface. In the low frequency region 
(\u — E a \ ~ \u> — fi\ -C min^ |A Q ^|), the diagonal ele- 
ments of the spectral functions for the other bands are 
expected to have negligible weight. This is because in 
equilibrium, in the non- interacting limit, A k ^i is sharply 



peaked at an energy separated from the chemical poten- 
tial by the direct band gap. With interactions, or out of 
equilibrium, there may be some incoherent weight around 
uj ~ /i but it is expected to be extremely small. Thus the 
other bands can also be eliminated, so that we obtain 
independent Keldysh equations for the diagonal spectral 
functions. 

Following this prescription, we obtain the differential 
equation C w {A a ) + C^(A a ) = with [H 

£ (1) (A») = (d Xj (L d!(X -M a ))(d Qj A a )-{X j ^Q j }, 

£W(A a ) = l{d Xj L d , a )(d Qj Q XkQh )A a 

-(d Xk L d , a p% kX .(d Qj A a ) - {Xj «- Qj} 
~{X k ^Q k } + {X k ,X ] ^Q k ,Q ] ). (8) 

Here the Berry curvatures associated with the a-th band, 
i.e. &>Q jXk and its 3 counter parts, are defined in terms 
of the renormalized gauge fields introduced in Eq. |7|. 
e.g. 

A4 a in Eq. ([5]) denotes an energy correction of the 1st or- 
der gradient expansion: M a = § J2j=o Hp^cXAx^ap ' 
(L d ,a — Ld.p) ■ [-4.qJ/3q + c.c. The Zeeman coupling en- 
ergy between orbital momentum and external magnetic 
field, which was previously found in the non-interacting 
case is encoded into the Y^=i terms in this energy 
correction. On the other hand, the temporal term, i.e. 
j = 0, describes a coupling of the electric dipole moment 
of quasi-particles to the external electric field, which is 
unique to an interacting FL [k|. 

We can now extract the low energy behavior of the 
spectral function. The renormalized energy dispersion of 
the a-th band e(q;X) is determined from 

L d , a -M a = e- E a (e, q; X) - M a (e, q; X) = 0. (10) 

Indeed, neglecting C^ 2 \ the form A a — S(L d ^ a — M a ), 
satisfies the Keldysh equation to the 2nd order in the 
gradient expansion. Therefore, without C^ 2 \ the spec- 
tral function is sharply peaked at lu — e, with a weight 
given by the conventional "renormalization factor" Z — 
[du}{L d ^ a — A^ Q )]|^, 1 =C . However, due to the presence of 

£( 2 ) , which contains terms proportional to A a rather 
than its derivative, the solution up to the 2nd order of 
the gradient expansion acquires an additional renormal- 
ization factor Z' = 1 — hC^x k Q k )|w=ei and 

A a = Z'8{L d>a - M a ) = ZZ'5(u - e). (11) 

To unmask the meaning of this new term Z' , let us 
specialize to the situation where the disequilibrium is 
generated by a physical (i.e. external) electromagnetic 
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gauge field (ao,a), with corresponding physical electro- 
magnetic fields b — V x a, e = V^ao — cVa. In this 
case, energy e(q,X) = eo(q + a (A)) + &o(-X"), where 
eo is the (renormalizcd) band energy in equilibrium. 
Then, introducing canonical momentum and frequency 
k = q+a(T, R) and J = u-&q(R) [Hj], Eq. $TT]) becomes 
A a = ZZ'S(u>' — eo(k)). Correspondingly, the Berry cur- 
vature appearing in the additional renormalization factor 
Z' are expressed in terms of the following artificial elec- 
tromagnetic fields (a.e/m.f.) estimated at u> — e, 



)\u 



(12) 



Namely, Z' simply reduces to an inner product be- 
tween real electromagnetic fields and these a.e/m.f, as 
in Eq. ((3]) . This result indicates a renormalization of the 
"angle resolved density of states" when both physical and 
a.e/m.f. are present. A similar modification (dependent 
upon B a ■ b only) was proposed in Ref. \\A . 

This theoretical observation leads us to propose the 
momentum-resolved detection of a.e/m.f. by ARPES, by 
observing the modification of the spectral weight by small 
applied electromagnetic fields. This is no doubt quite 
difficult in practice, but demonstrates that the systematic 
measurement of the artificial fields is possible in principle. 

We now turn from the spectral function to the quasi- 
particle dynamics, by deriving the effective Boltzmann 
equation for quasi-particlcs in a U{1) FL. We start with 
the observation that the dissipationless Keldysh equa- 
tion in Eq. §5§ holds also for the lesser Green's func- 
tion, g < (l,l / ) = i(ip a ,(r' ,t') ■ ip a (r,t)), under the same 
assumptions [lfj |. Applying the same unitary transfor- 
mation as above, i.e. g K = W ■ g < • U, one then ob- 
tains exactly the same reduced Keldysh equation for the 
(a, a)-component of g K as in Eq. ([5|). This lesser Green's 
function may be decomposed according to g<(Q;X) = 
iA a (Q;X) ■ f a (Q; X) where f a (Q;X) is a generalized 
Fermi distribution [10]. Inserting this into the reduced 
Keldysh equation thus obtained, we can readily obtain 
the EOM for /„, 

= ZZ'S(lu - e) x {d X] (L d . a - M a )d Q] f a 

- dx k n% hXj d Qj f a + d Qk n XkXj d Qj f a - (Xj r<-> Qj)}. 

Because of the sharp peak in A a ~ 5(u> — e), when inte- 
grating this EOM w.r.t. ui, we get an effective Boltzmann 
equation for a quasi-particle occupation number in q — R 
space; n a (q : R,T) = f a (q,e : R,T), 



{d q 



drn a 



TRj 



(l-Q^ + ^.e) n« e 

= (d Rj e+(d T e) 

+ (d qk c) V Q RkR] ) 3 qi n a - { Qj ^R,}. (13) 

Note that, although all the curvature terms in the above 
are now "on-shell" quantities, the partial q, R and T- 
derivatives encoded there applies only on their explicit 
dependence and do not apply on their arguments of e. 



In the presence of external electromagnetic fields, this 
effective Boltzmann equation can be substantially simpli- 
fied in terms of canonical momentum introduced above, 



= drn a + — e + b x v 
|bx(exB^)-bx((bxv)x B'J 



■ d k n 



v + e x B' - (b x v) x B' 



■ d R n a: 



(14) 



where we used abbreviated notations B' a = B a — £ a x v 
and v = c^e fljj . The partial X-derivative of n a was 
taken with the canonical momentum k fixed. By equat- 
ing this with the continuity equation in phase space, 
= dTn a + (drk) ■ d^n a + (drR) ■ d R n a , we obtain 
the quasi-particle equation of motion, Eq. ([1]), to the 
accuracy of the 2nd order gradient expansion. Further- 
more, extracting the conserved current from the conti- 
nuity equation in zero external field (b = 0), one finds 
j a = f k [\r + e x B' a ]n a . The usual solution of the Boltz- 
mann equation in linear response therefore gives imme- 
diately the anomalous Hall conductivity in Eq. |(2J). 
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